Nanometer-scale spatial modulation of an inter-atomic interaction in a Bose-Einstein 

condensate 
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We demonstrate nanometer-scale spatial control of inter-atomic interactions in a Bose-Einstein 
condensate of ytterbium(Yb). A pulsed optical standing wave, tuned near an optical Feshbach 
resonance varies the s-wave scattering length continuously across the standing wave pattern. The 
modulated mean-field energy with a spatial period of every 278 nm is monitored by a diffraction 
pattern in a time-of-flight image. We observe a wide scattering length control of up to 160 nm. The 
demonstrated spatial modulation of the scattering length proves that the high resolution control of 
atomic interactions is possible. 
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Tunability of interatomic interactions using Feshbach 
resonance methods [l( has opened up a variety of possi- 
bilities in the study of ultracold quantum gases such as 
formation of ultracold molecules [2j , a Bose-Einstein con- 
densate (BEC) to Bardeen-Cooper-Schrieffer crossover 
with fermionic gases [3| , simulation of super- nova (Bosc- 
nova) [J], and Efimov trimer states [5]. In these stud- 
ies, the Feshbach resonances are induced by magnetically 
tuning the hyperfine energy level of a bound state to the 
vicinity of the incoming scattering state. The length scale 
for application of the Feshbach field has, so far, been 
much larger than the size of the atomic sample. 

Many interesting experiments such as a sonic-analogue 
of black holes IfiJ, quantum simulation of novel Hub- 
bard model [7|, |8| , and unusual Bose-nova phenomena [9( 
could be explored with spatial modulation of the inter- 
atomic interaction on short length scales. An alternative 
approach, an optical Feshbach resonance (OFR), utiliz- 
ing optical coupling between bound and scattering states 
10], is a promising technique for providing this possibil- 
ity. With OFR methods, the intensity and the detuning 
of the coupling laser are used to control the scattering 
length llj. Fine modulation of laser intensity using opti- 
cal standing wave and holographic technique is routinely 
shown in other applications [12l . |13| , and these techniques 
combined with OFR should realize the fine spatial mod- 
ulation of the scattering length. Inherent nature of the 
optical control also enables fast manipulation and, more- 
over, optical transitions can be found in almost any states 
and the control is possible for numerous situations and 
atomic species. 

Important progress along this direction has been re- 
ported in Ref. [14j where successful control of a mag- 
netic Feshbach resonance of alkali-metal atoms with a 
laser light was demonstrated. The fine spatial control of 
the scattering length has, however, not yet been demon- 
strated, so far. In this Letter, we report our successful 
demonstration of OFR modulation of the s-wave scatter- 
ing lengths on the scale of hundreds of nm. 

For the demonstration of submicron control of the scat- 



tering length, we apply a pulsed optical standing wave to 
a BEC of ytterbium ( m Yb) atoms. The pulsed light for 
OFR is tuned to the vicinity of the weak intercombina- 
tion 1 So- 3 Pi photoassociation (PA) resonances with the 
resonant wavelength of 556 nm, which results in a modu- 
lation of the mean-field energy of the condensate with a 
spatial period of 278 nm. The resulting phase modulation 
of the condensate is observed in the diffraction pattern 
in a time-of-flight (TOF) images. The observed behav- 
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FIG. 1: (color online), (a) Schematic of the experimental 
setup. A 174 Yb condensate is irradiated with the standing 
wave formed by an OFR laser. The diffraction pattern in 
the time-of-flight (TOF) image is observed, (b) Energy dia- 
gram of the relevant states for the experiment. The 1 So- 3 Pi 
photoassociation transitions to the vibrational states v' = 11, 
12, and 13 are used for the OFR. (c) Typical diffraction pat- 
tern obtained in the experiment. Each peak in the image 
represents different momentum components imparted by the 
pulsed lattice, (d) An integrated column density of the image 
where the red dots and a blue line corresponds to the data 
and the fitted line for the determination of ft, respectively. 



ior is well explained by the semi-classical theory of Bohn 
and Julienne |15| . Our analysis shows scattering length 
variation of up to 160 nm over a comparable distance of 
only 280nm. 

To determine the expected change of scattering length 
a and PA rate K due to the applied laser field, the for- 
malism of Bohn and Julienne is valid under our condi- 
tions. In their formalism, both parameters can be cal- 
culated in terms of the radiative decay rate of molecule 
7, the light-induced width T and the detuning of the 
OFR laser from a photoassociation resonance A. The 
energy level diagrams relevant to this work is shown 
in Fig. 1(b). The light-induced width T is given by 
T(I) = 3IjX 3 frotfFc/8nc, where f rot = 1/3 is the rota- 
tional factor for the relevant transition, fpc is the Frank- 
Condon factor, c is the speed of light in a vacuum, and 
A = 555.8 nm and / are the wavelength and the intensity 
of the OFR laser, respectively. For the weak excitation 
regime, where r <C 7, the scattering parameters are given 
by 
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where a bg = 5.55 nm is the background s-wave scattering 
length of 174 Yb, /i = to/2 is the reduced mass of two 

scattering atoms with mass to, and k — '' J ' ' 



2R TF ^ 

the wave number calculated from the momentum of the 
condensate with Thomas- Fermi radius Rtf |16| - 5a is 
the variation of the scattering length due to the OFR 
and is proportional to /. 

Despite its potentially wide applicability, the use of 



OFR has been scant (ll|. A reason for this is that usu- 



ally the optical coupling also induces inelastic scattering, 
leading to rapid atom depletion. A use of narrow transi- 
tions in alkaline-earth-metal-like atoms 17| to avoid in- 
elastic scattering loss was successfully demonstrated in 
a recent experiment in our group using the ^o^Pi in- 
tercombination transition in thermal gases of 172 Yb and 

176 Yb - 
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In that earlier work a relatively low inelastic scatter- 
ing rate was observed, an order of magnitude less than 
the case of alkaline atoms indicating the possibility of 
observing large tunings of scattering length over short 
distances, as we report here. In this work we extend this 
technique to a condensate of 174 Yb. 

To demonstrate the capability of OFR to modulate 
the s-wave scattering length on short length scale and 
also on short time scale, we use a pulsed optical lattice 
beam as the OFR light, as shown in Fig. 1(a). The appli- 
cation of a pulsed optical lattice beam generally results 
in a diffraction of a released BEC, which is previously 
studied both theoretically and experimentally |12]. In a 



phase modulation regime, where the exposure time of the 
optical lattice r is much smaller than the minimum clas- 
sical oscillation period of the formed lattice, the effect of 
the lattice is treated as a thin grating. The phase of the 
condensate modified by the lattice has a form 



(*) = 



U(z)t 



(3) 



with U(z) = Uocos 2 (qz), where q is the wavenumber of 
the lattice laser and z is the direction along the lattice 
laser propagation. The potential J7o experienced by the 
atoms in the ground state at an anti-node of the optical 
lattice is given by 
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The first term describes the atomic light shift, where 
£1 and A at om corresponds to the Rabi frequency and the 
detuning of the OFR laser to the atomic excited state, re- 
spectively. The second term represents the OFR-induced 
shift of the mean-field energy, and thus is proportional 
to the scattering length variation Sa and the atom den- 
sity n a . It is noted that there is a mean-field energy 

4irH 2 n a 
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L ab g across the condensate. However, it is 



not susceptible to the OFR laser and does not contribute 
to the diffraction pattern. It is clear from Eq.([T]) and 
flU), the OFR dispersively varies the scattering length, 
and therefore the mean-field energy, across a photoasso- 
ciation resonance. From the diffraction pattern of the 
condensate generated from the phase modulation tfi(z) 
imparted by the OFR, one should be able to extract the 
variation of the scattering length. 

The method for the all-optical formation of 174 Yb con- 
densate is described in Ref. [19]. After the evaporation 
in a crossed far-off resonant trap (FORT), an almost 
pure 174 Yb condensate is prepared with an atom num- 
ber of up to 1.5 xlO 5 in a typical harmonic trap potential 
Gj = (ojxUjyUJz) 1 / 3 = 2n x 130 Hz. A simple schematic 
of the experimental process after the preparation of con- 
densate is shown in Fig.flja). Following the condensate 
formation, we release the condensate from the trap by 
turning off the FORT lasers. At the release, the OFR 
laser pulse of a ID optical lattice is turned on for a sev- 
eral microseconds with a typical power of 1-100 fiW. The 
OFR laser is tuned near the 1 So- 3 Pi photoassociation 
resonances with the vibrational quantum number v' = 
11, 12 and 13, which correspond to the detuning A atom 
of -69, -117, -192 MHz, respectively 0. After a TOF- 
time of typically 10 ms, the absorption image is taken for 
the diffraction pattern analysis. For every OFR pulse, 
the power of the pulse is monitored by a fast photodiodc 
and recorded with an oscilloscope to compensate for the 
power instability. 

A typical image and a column density of the obtained 
diffraction pattern are shown in Fig.QJc) andQJd), re- 
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FIG. 2: (color online), (a) Observed modulation index ft, 
(b) number of atoms remained after the OFR pulse, and (c) 
calculated scattering length for the v' — 12 transition. We 
first fit the variation of j3 using the formalism by Bohn and 
Julienne to determine light-induced width T. Blue lines in 
(b) and (c) are calculated atom loss and scattering length 
variation using the F obtained from the fit in (a), respectively. 



spectively. Each peak in the image represents the mo- 
mentum component p n — n2hq, (n = 0, ±1, ±2, . . .) im- 
parted by the pulsed lattice. Population distribution for 
each momentum state is described as P n — J%(/3), where 
j3 = Uqt/2H and J n (z) are Bessel functions of the first 
kind [12J. For the image analysis, we first sum the image 
along the vertical axis to obtain an integrated column 
density profile in the horizontal axis, and fit the density 
profile by the Bessel functions to obtain the modulation 
index /3. Each peak is convoluted with a Gaussian func- 
tion to account for the momentum spread and the finite 
image resolution. It is also noted that we subtract the 
background distribution in the images. There are two 
origins of the background distribution; the residual ther- 
mal component and the elastic scattering of wave packets 
of condensate with different momentum during the TOF 
22j. Strictly speaking, the atom distribution from the 
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later phenomena are non-Gaussian. However, due to low 
atom counts from the elastic scattering we do not observe 
clear deviation from the Gaussian fit for the background 
distribution. 

The images of the diffraction pattern across the pho- 
toassociation resonances are collected by scanning the 
OFR laser frequency. We keep the OFR pulse length to 
be short and all the diffraction pattern images collected 
are in the phase modulation regime. Shown in Fig. [5] 
are typical /3 and atom number after the OFR pulse as 
obtained from the images. The observed /? has a disper- 
sive shape centered at the photoassociation resonance, 
where atom loss is maximum. The slight increase in j3 
towards the atomic resonance is due to the variation of 
the atomic light shift. To extract the scattering length 
from the data, we first fit the observed /3 to determine 
r and f2. We then calculate the atom loss and the vari- 
ation of the scattering length 5 a from these parameters. 
Atom loss obtained from the data agrees well with that 
obtained from the fit to the /3 variation. 

The clear dispersive variation in (3 is direct evidence of 
the OFR effect. We further investigate the behavior of 



FIG. 3: (color online), (a) Power dependence of the scatter- 
ing length variation along with a linear fit line. The linear 
dependence is consistent with the weak excitation expected 
from the V value obtained from the data, (b) Atom num- 
ber dependence of 5/3 along with the best fitted curve in blue. 
The OFR induced mean-field variation is expected to depends 
on the atom density which is proportional to N 2 ' 5 in a har- 
monic potential. From the fit, we obtained p — 0.41 for N p 
dependence fitting, which agrees well with the predicted den- 
sity dependence. The error bars in both figures show the 
root mean squared error of the data obtained from the fitting 
procedure of the modulation index /3. 



the signals to confirm the OFR effect. First, the variation 
of Sa with respect to the OFR laser power is measured. 
The fact that r/27r < 20 kHz for all the data obtained 
means that all the measurements are performed in the 
weak excitation condition T ^; 7 (=27rx364 kHz) where 
5a should be proportional to the laser intensity as shown 
in Eq. (1). The results shown in Fig.^a) clearly confirm 
the linear dependence, which is consistent with the weak 
excitation regime. The observed variation of the scatter- 
ing length up to 160 nm is quite large. It is important 
to note that 5a varies continuously from the anti-node, 
where the effect of OFR is maximum, to the node where 
the effect of OFR is absent. The OFR laser wavelength 
of A = 555.8 nm results in a periodic modulation of 5a 
in approximately 278 nm spacing. 

Second, the variation of the modulation index across 
the photoassociation resonance 5f3 with respect to the 
atom number is measured. The OFR effect is propor- 
tional to the density as is shown in Eq. @. For a given 
trap frequency, the density n a is proportional to TV 2 ' 5 
23] . Figure[3Jb) shows the variation of 5/3 as a function 
of the atom number with a best fit curve of N p shown 
with a blue line. From the best fit, we obtain p = 0.41 
and the value agrees well with the 7V 2 / 5 -dependence pre- 
dicted. The agreement between the data and the curve 
clearly shows that the observed 5(3 comes from the mean- 
field energy of the condensate. 

The effect of the OFR depends on the Frank-Condon 
factor of the photoassociation resonance chosen. From 
the same measurement on different photoassociation res- 
onances, we obtain the optical length defined as L opt = 
r/27 which describes the strength of the OFR coupling 
for each resonance. We plot the optical length obtained 
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FIG. 4: (color online). Comparison between calculated (red 
dots) and measured (blue dots) L opt . The rapid increase of 
L op t for the shallower detuning are observed for both results. 
The discrepancy between two results may be due to the un- 
certainty in the laser intensity estimation at the condensate. 



from the measurement (blue dots) and the numerical cal- 
culation (red dots) |24j in Fig.0J As the transition gets 
shallower toward the atomic resonance, the measured op- 
tical length shows faster increase in the optical length 
than the numerically calculated values. The tendency in 
the variation of L opt , however, is consistent between two 
results. Part of the discrepancy may be due to the uncer- 
tainty in the laser intensity estimation at the condensate. 

In conclusion, we demonstrate rapid, nanometer-scale 
modulation of the scattering length using optical Fesh- 
bach resonances in a m Yb condensate. A large scatter- 
ing length variation of up to 160 nm is obtained with a 
spatial modulation of 278 nm by observing a diffraction 
pattern in a time-of-flight image. The demonstrated fine 
spatial modulation of the scattering length opens up a 
wide variety of applications. New applications may ben- 
efit not only from the fine spatial resolution offered by a 
OFR, but also the inherently available fast manipulation 
[25| . The time constant for the manipulation of the OFR 
can be quite fast compared to the magnetic Feshbach 
resonances and could be used to control and study the 
non-equilibrium dynamics of the condensate. In addi- 
tion to these possibilities, creation of multiple frequency 
super-lattices with frequency components tuned to dif- 
ferent OFR resonances may create intricate patterns of 
large variation in the interaction strength. Independent 
control of different scattering lengths [26J] , and local ma- 
nipulation of the condensate or the condensate in the 
lattices may also be possible. 
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